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摘 要 : 图 像 模拟 在 天 文学 研究 中 已 经 发 挥 越 来 越 重要 的 作用 ， 其 重要 性 主要 
表现 在 以 下 几 个 方面 : 通过 模拟 可 以 为 准备 研制 的 天 文 观测 设备 评估 提供 重要 的 
依据 ; 通过 对 模拟 数据 的 处 理 可 以 验证 数据 处 理 的 算法 。 本 文 对 图 像 模拟 的 整个 
过 程 进 行 描述 ， 并 在 计算 上 对 图 像 模拟 的 局 部 方法 进行 优化 ， 使 模拟 结果 能 够 高 
效 准 确 的 生成 。 本 文 所 做 的 主要 工作 包含 以 下 几 个 方面 : 对 天 体 的 形态 的 模拟 ， 
包括 点 源 (如 恒星 ) MRE (星系 ) ， 而 这 里 主要 是 对 展 源 建 模 ; 模拟 观测 条 件 ， 
主要 是 仪器 和 天 光 背 景 等 产生 的 噪声 ; 模拟 点 扩散 画 数 (PSF) ; 从 计算 和 程序 上 
的 方法 对 计算 流程 进行 优化 ; 对 结果 进行 分 析 。 本 文中 所 条 用 的 数据 是 哈 勃 超 深 
空 场 (HUDF) ACS WFC i (F775) 波 段 的 数据 ， 其 极限 星 等 达到 
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天 文 图 像 模拟 是 将 图 像 处 理 的 技术 与 天 文学 模型 相 结 合 ， 这 种 模拟 在 天 文学 研 
究 中 已 经 越 来 越 广泛 的 使 用 。 模 拟 有 遇见 未 来 之 功效 ， 通 过 模拟 可 以 对 未 来 的 设 
备 进行 更 加 科学 的 评估 。 基 于 这 种 原因 ， 越 来 越 多 人 开始 从 事 天 文 图 像 模 拟 方面 
的 研究 。 在 对 天 文 的 观测 研究 中 ， 有 一 些 被 学 者 和 观测 者 熟识 的 天 文 观测 模拟 软 
件 ， 如 WorldWide Telescope, Astronomy CCD Calculator, Stellarium, X £E 
软件 对 天 体 空间 位 置 、 天 体 特性 描述 的 非常 细致 ， 给 观测 者 带 来 了 很 大 的 便利 ， 
但 是 对 于 遥远 星系 的 描述 这 些 软件 还 不 是 很 完善 。 随 着 大 型 观测 设备 的 出 现 ， 人 
们 对 未 知 天 空 的 观测 变 的 越 来 越 琐 远 ， 星 系 模拟 也 成 为 评估 大 型 天 文 设备 很 重要 
的 一 部 分 。BertinD 在 其 文中 介绍 了 SkyMaker， 该 软件 很 细致 的 描述 了 点 扩散 函 
数 ， 包 括 了 大 和 气 扰动 、 望 远 镜 运 动 、 像 差 、 衍 射 等 各 个 因素 ， 对 于 星系 的 形态 其 
利用 了 De Vaucouleurs 形态 。Refregier “等 人 利用 shapelet 的 方法 对 星系 的 形态 
进行 模拟 ， 该 方法 对 星系 形态 的 模拟 采用 多 维 的 模拟 方法 ， 文 中 与 wavelet 等 方法 
进行 比较 ， 对 残 差 进行 分 析 ， 证 明 其 方法 的 可 靠 性 。LSST2 成 立 了 专门 的 图 像 模拟 
团队 ， 该 团队 主要 的 工作 就 是 为 LSST 做 模拟 工作 。 他 们 的 模拟 工作 是 针对 全 天 区 
的 模拟 ， 利 用 宇宙 模型 建立 星 表 来 确定 天 体 的 分 布 ， 并 且 建 立 了 多 种 天 体现 象 的 
模型 〈 天 体 模型 包括 恒星 、 星 系 、 小 行星 等 ) ， 在 传输 过 程 中 考虑 天 气 情况 、 大 
气 扰动 等 ， 最 后 还 考虑 了 结合 巡天 设备 影响 所 产生 的 PSF 函数 。 利 用 准确 的 模拟 结 
果 可 以 对 LSST 系统 的 各 项 性 能 进行 分 析 处 理 。 另 外 ，LSST 对 天 体形 态 的 模拟 是 采 
用 了 蒙特 卡 罗 的 方法 ， 模 拟 光子 的 射 入 。Joel Berge 等 发 表 的 文章 中 提 到 了 模 
拟 的 两 种 不 同方 法 ， 一 种 是 类 似 LSST 模拟 光子 射 入 的 方法 ， 另 一 种 是 按照 像素 模 
拟 的 方法 ， 并 比较 两 种 方法 的 差异 ， 该 文章 用 到 了 简单 的 数学 模型 ， 主 要 强调 了 
处 理 的 迅速 。 
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本 文中 所 做 的 模拟 是 利用 Hubble 望远镜 的 数据 和 观测 条 件 。 文 中 的 图 像 模 拟 
所 做 的 工作 包含 : 1) 模拟 天 体 分 布 和 天 体 的 形态 ， 主 要 是 指 星系 的 形态 2) 模 
拟 观测 条 件 以 及 各 种 噪声 ， 包 括 天 光 背 景 、 暗流 、 读 出 噪声 等 ，3) 模拟 仪器 的 
PSF; 4) 对 观测 结果 进行 统计 分 析 。 图 1 为 本 文 图 像 模 拟 的 流程 图 。 
本 文 的 主要 目标 是 利用 简单 通用 的 模型 对 天 文 图 像 实现 快速 的 模拟 。 在 程序 编写 上 应 用 更 
多 的 优化 ， 使 程序 达到 快速 可 用 的 目的 。 本 文 还 将 利用 SExtractor 对 模拟 结果 进行 提取 ， 将 
提取 的 结果 与 参考 结果 进行 比较 分 析 ， 通 过 分 析 确 定 模 拟 结 果 在 实际 应 用 上 的 可 行 性 。 本 文 1 
部 分 介绍 图 像 模拟 的 几 个 关键 模型 和 方法 ， 主 要 包括 星系 形态 的 模拟 模型 、PSF 模型 以 及 观测 
条 件 模 拟 等 ，2 部 分 是 对 模拟 结果 的 分 析 ; 3 部 分 对 整个 方法 流程 总 结 。 


图 1 图 像 模拟 流程 图 


Fig. 1 Flow Diagram of Image Simulation 


1 图 像 模拟 的 关键 模型 及 方法 


本 文中 的 模拟 所 利用 的 数据 是 HUDF ACS WFC i (F775) 波段 的 数据 ， 该 数据 生 
成 的 图 像 约 为 10000*10000 个 像素 ， 由 于 图 像 在 空中 存在 一 定 的 旋转 ， 图 像 四 个 角 
包含 未 拍摄 的 区 域 ， 本 文 为 了 实验 方便 ， 裁 取 了 其 中 4500*4500 个 像素 的 中 间 区 域 ， 
面积 约 为 5.1 arcmin*。 文 中 的 模型 采用 了 天 文学 中 一 些 简单 高 效 的 模型 ， 主 要 目的 
就 是 为 了 对 数据 进行 快速 处 理 。 在 观测 条 件 上 ， 模 拟 了 Hubble 望远镜 在 600KM 的 
高 空中 的 观测 条 件 。 


1.1 星 系 的 模型 


文中 的 星系 的 模型 采用 了 Sersic 模型 5, MERRITT: 
(1) 
式 (1) 中 表示 R 处 的 流量 密度 ， 表 示 有 效 半径 处 的 流量 密度 ， 为 有 效 半径 ，n 为 
Sersic 系数 ， 通 过 该 系数 可 以 控制 星系 的 弥散 程度 ， 为 常数 。 
通过 变换 可 以 将 (1) 式 子 变换 如 下 表示 方法 : 
(2) 
(2) 式 中 表示 中 心 位 置 的 流量 密度 ，n 为 Sersic 系数 ，a 为 常数 ， 表 示 当 R 为 a 时 
的 流量 密度 为 中 心 位 置 流 量 密度 自 


rea 


若 星系 为 圆 形 区 域 ， 对 (2) 式 在 圆 面 内 积分 ， 得 到 星系 的 流量 
(3) 
(3) 式 中 为 不 完全 gamma 函数 ，。 当 积分 的 半径 趋 于 无 穷 时 
(4) 
(4) 式 中 为 gamma 函数 。 
在 模拟 中 ， 通 常用 椭 率 来 表示 星系 的 形状 ， 这 就 需要 将 星系 的 流量 密度 分 布 在 


整个 椭圆 平面 内 ， 定 义 椭圆 面 的 椭 率 为 e， 则 (2) 式 中 R 可 以 用 归 一 化 的 椭圆 半径 表 
人 小 
(5) 
(5) 式 中 表示 坐标 为 ， 坐 标的 原点 为 椭圆 的 中 心 。 
星系 在 图 像 中 的 表现 ， 除 了 弥散 程度 和 形状 的 特性 外 ， 还 存在 着 方向 性 ， 这 就 
需要 对 已 经 计算 出 来 的 椭圆 形 的 Sersic 分 布 进行 坐标 变换 ， 假 设 星系 坐标 系 〈 长 轴 
为 x 轴 ， 短 轴 为 y 轴 ) 与 图 像 坐标 系 的 夹 角 为 ， 则 变换 矩阵 为 


则 新 的 坐标 位 置 为 
(6) 

通过 (1)-(6) 的 计算 可 以 模拟 出 星系 的 形状 、 大 小 、 方 向 等 要 素 ， 而 对 星系 的 模 
拟 所 需要 的 人 参数， 如 (Sersic 系数 ) 、 (有 效 半径 ) 、 亮 度 等 参数 都 是 通过 对 
HUDF ACS WFC 的 数据 进行 提取 获得 。 在 Coe D. 等 四 的 文中 利用 SExtractor 对 HUDF 
的 数据 进行 提取 ， 并 和 Beckwith 提供 的 数据 进行 对 比 ， 结 果 比 较 理 想 ， 文 中 所 有 
用 到 的 参数 均 来 自 于 Coe D. 等 所 提取 的 星 表 。 
1. 2 观测 条 件 及 噪声 的 模拟 

对 于 观测 条 件 的 模拟 ， 在 图 像 中 表现 其 实 就 是 对 噪声 的 还 原 ， 还 原 的 真实 与 否 
是 影响 模拟 图 像 效 果 的 最 主要 因素 。 在 观测 中 所 产生 的 噪声 主要 包括 固有 了 噪声 、 
言 号 噪声 。 

固有 了 噪声 是 指 源 于 探测 器 自身 的 噪声 ， 包 括 热 噪声 、 散 粒 噪声 、 生 成 -复合 噪 
声 和 内 变 噪声 ， 在 文中 用 读 出 噪声 将 其 表示 。 为 了 和 Hubble 的 数据 进行 对 比 ， 读 
出 噪声 是 利用 Hubble 的 曝光 时 间 计 算 器 外 计算 得 出 。 

信号 噪声 在 文中 主要 包含 的 是 背景 噪声 和 由 于 信号 的 量子 性 所 产生 的 随机 噪声 。 
背景 噪声 利用 了 HST 测 得 的 天 光 背 景 噪声 ， 包 括 黄道 光 和 赤道 光 。 

对 于 信和 号， 由 于 其 量子 性 ， 到 达 探 测 器 会 有 一 定 的 随机 性 ， 基 于 这 种 特性 也 会 
产生 一 定 的 噪声 ， 在 文中 用 泊 松 分 布 表示 这 种 特性 。 
1. 3 点 扩散 函数 模拟 

对 于 地 面 望远镜 ， 影 响 点 扩散 函数 的 主要 因素 有 大 气 扰动 和 仪器 自身 的 光学 特 
性 。 而 对 于 Hubble 望远镜 这 样 的 空间 设备 来 说 ， 观 测 条 件 要 远 远 好 于 地 面 望 远 镜 ， 
因此 对 于 点 扩散 函数 的 影响 仅仅 是 来 自 观测 设备 的 光学 特性 。 

文中 所 做 的 模拟 是 对 空间 望远镜 成 像 模 拟 ， 点 扩散 函数 采用 Gaussian 函数 表示 ， 
如 下 : 

(7) 
上 式 中 ，16 为 常数 ， 表 示 中 心 位 置 的 流量 密度 ， 可 以 控制 Gaussian 函数 的 幅度 ， 
[x, y AERE, o 用 来 控制 图 像 的 质量 ， 对 于 Gaussian 函数 通常 可 以 通过 所 
给 的 半 高 全 宽 (FWHM) 所 求 得 。 
(8) 

而 对 于 地 面 的 望远镜 通常 需要 考虑 大 气 扰动 的 效应 ， 所 以 模拟 中 所 用 到 的 点 扩 

散 函 数 需要 考虑 观测 设备 和 大 气 扰 动 的 综合 效应 ， 对 于 这 类 点 扩散 函数 模拟 通常 


采用 Moffat RA, 如 下 : 
(9) 
其 中 ，lo 表 示 中 心 位 置 处 的 流量 密度 ，a M8 为 尺度 因子 ， 可 以 控制 函数 的 弥散 
程度 ，[x, y] 为 图 像 坐 标 ，a MB 关系 如 下 : 
(10) 
文 [10] 中 提 到 ， 当 8 =4. 765 时 ，Moffat PSF 可 以 很 好 的 模拟 大 气 扰动 。 


1.4 子 像素 级 处 理 


模拟 的 其 中 一 个 目的 就 是 对 未 来 某 一 个 观测 设备 的 评估 ， 文 中 所 模拟 所 利用 的 
数据 只 符合 Hubble ACS 的 特性 ， 例 如 像素 的 大 小 为 0.03X0.03 
arcsec (Drizzling 之 后 的 像素 大 小 ) ， 为 了 使 程序 能 够 为 其 他 观测 设备 评估 所 利用 ， 
文中 的 模拟 需要 做 到 子 像素 级 。 至 于 所 要 做 到 的 子 像素 级 别 需 要 根据 具体 的 情况 
而 定 。 

本 文 所 做 的 子 像素 级 处 理 主要 是 应 用 在 天 体形 态 、 亮 度 的 模拟 过 程 中 ， 根 据 需 
要 将 原 像 素 分 成 若干 个 子 像素 若干 个 子 像素 能 将 原 像素 均 分 ) ， 根 据 天 体 的 特 
性 ， 在 子 像素 中 计算 得 到 该 像素 的 值 。 这 一 计算 ， 单 一 看 来 十 分 简单 ， 但 综合 起 
来 ， 计 算 量 会 变 得 十 分 庞大 ， 这 就 需要 对 程序 进行 其 他 优化 。 


1.5 并行 计算 


eG xs、 亮 度 的 要 素 时 ， 由 于 将 原 有 像素 拆 分 为 子 像素 计算 ， 这 样 

导致 计算 量 倍增 。 另 外 ， 每 一 个 天 体 都 是 独立 的 个 体 ， 单 独 计算 对 其 他 天 体 并 

会 产生 任何 影响 。 所 以 ， 文中 将 计算 每 二 个 天 体 的 形态 、 亮 度 等 要 素 进行 独立 
计算 这 样 ， 将 模拟 中 的 所 有 天 体 按照 CPU 的 数目 进行 平均 分 配 ， 不 同 的 CPU 计算 
完 之 后 的 结果 返 给 主 进程 进行 合并 ， 形 成 最 终 模拟 的 整 幅 图 像 。 


2 结果 分 析 


2. 1 参考 样本 分 析 


文中 用 到 Coe D. 等 提供 的 数据 ， 该 数据 中 包含 了 提取 出 的 Sersic 系数 ， 而 
Beckwith 提供 了 一 份 更 加 准确 的 数据 ， 不 过 在 该 数据 中 对 于 天 体 弥散 程度 的 描述 
只 给 出 了 Gaussian 的 参数 ， 但 是 其 位 置 和 亮度 要 更 加 准确 ， 事 实 上 ，Coe D. 等 提 
供 的 数据 中 也 给 出 了 一 些 与 Bechwith 的 数据 的 对 比 ， 例 如 ， 前 者 的 数据 中 给 出 了 
中 心 位 置 与 Bechwith 的 对 比 ， 本 文 对 两 个 数据 做 了 更 详细 的 对 比 。 

图 2 是 对 两 份 数据 提取 天 体 的 亮度 的 对 比 ， 蓝 色 的 线 为 Coe D. 等 的 数据 天 体 
亮度 的 统计 ， 红 色 的 线 为 Bechwith 的 数据 天 体 亮度 的 统计 。 从 该 图 中 可 以 看 出 ， 
两 份 数 据 的 亮度 存在 着 一 定 偏差 ，Coe D. 等 提供 的 数据 天 体 亮度 要 比 实际 偏 大 。 
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图 2 Coe D. 等 提供 的 数据 与 Bechwith 的 数据 天 体 亮度 对 比 


Fig. 2 AB magnitude comparison between Coe et al. catalog and Bechwith et al. catalog 


图 3 是 对 两 个 样本 有 效 半径 的 对 比 结果 ， 横 坐标 是 提取 天 体 的 亮度 信息 ， 纵 坐 
标 是 Benchwi ch 的 数据 中 天 体 的 有 效 半径 与 Coe D. 等 的 数据 中 的 相应 的 天 体 的 有 
效 半径 之 差 。 从 该 图 中 可 以 看 出 ， 差 值 绝 大 部 分 都 集中 在 0 值 附近 ， 但 是 Coe D. 
等 的 数据 有 部 分 半径 会 偏 大 。 

综合 上 述 对 样本 的 对 比 ， 我 们 利用 Coe D. 等 人 提供 的 数据 所 模拟 的 图 像 与 真实 
图 像 会 存在 亮度 偏差 ， 另外， 对 于 天 体 的 大 小 也 会 有 部 分 天 体 与 真实 图 像 存 在 一 
定 差别 。 从 这 一 点 ， 也 可 以 解释 后 面 章节 中 实验 结果 图 4 的 模拟 结果 存在 一 些 比 
原 图 中 大 的 天 体 的 现象 。 
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图 3 两 种 参考 样本 提取 天 体 有 效 半 径 的 对 比 
Fig. 3 Effective Radius comparison between Coe et al. catalog and Bechwith et al. catalog 
2. 2 模拟 结果 及 分 析 


图 4 展示 了 模拟 的 结果 ，a 图 是 模拟 图 像 ，b 图 是 drizzling 之 后 原始 图 像 。 
由 于 原 图 在 空间 中 存在 一 定 角 度 ， 并 且 四 个 角 上 没有 目标 信息 ， 所 以 本 文 所 用 到 


的 数据 和 图 像 是 对 原 图 的 和 截取， 截取 了 原 图 x 方向 3001 到 7500 像素 ，y 方向 3001 
到 7500 像素 的 方形 区 域 ， 总 共 4500*4500 个 像素 ，drizzling 之 后 每 个 像素 的 视 
场 约 为 0. 03X0.03 arcsec-， 模 拟 区 域 的 视 场 为 5.1 arcmin’. HST 团队 给 出 ACS 
WFC i (F775) 波段 零点 星 等 (zeropoint) 的 值 为 25. 654， 这 里 零点 的 意思 就 是 
流量 为 1 e- 时 天 体 的 亮度 值 。 根 据 式 子 (11) 可 以 计算 出 天 体 的 流量 
(11) 

由 式 (11) 所 计算 出 的 流量 值 是 信号 的 流量 。 而 噪声 包括 天 光 背 景 、 暗 流 、 读 出 噪 
声 等 ， 是 利用 hubble 提供 的 曝光 时 间 计 算 器 计算 得 到 。 计 算 的 结果 为 以 e- AB 
位 的 流量 。 所 有 计算 的 结果 除 以 增益 〈 单 位 ，e-/DN) 就 得 到 表现 在 图 像 上 的 数字 
值 。 

HUDF ACS WFC | 波段 的 极限 星 等 为 29 等 ， 所 以 本 文 提 取 的 截止 亮度 到 29 等 。 
文中 用 SExtractor 作为 提取 工具 ，SExtractor 的 提取 主要 参数 见 附件 1。 表 1 是 
对 提取 数目 的 统计 ， 人 参考 的 对 象 为 Coe D., Benitez N. 等 提取 的 数据 。 表 中 有 四 
项 数据 ， 第 一 项 是 对 原 图 提取 数据 的 统计 ， 第 二 项 是 对 模拟 图 像 提取 数目 的 统计 , 
第 三 项 是 粗略 的 做 了 匹配 统计 ， 匹 配 的 方法 是 将 模拟 数据 提取 出 来 的 天 体 的 中 心 与 
Coe D., Benitez N. 等 提取 的 数据 的 天 体 的 中 心 作 比较 ， 选 出 距离 最 接近 的 一 对 
被 认为 是 同一 个 天 体 。 第 四 项 是 计算 匹配 率 的 比例 ， 该 比例 通过 式 (12) 计算 得 到 

(12) 
上 式 中 ref num 是 参考 数据 提取 天 体 的 数目 ，matched sim num 是 对 模拟 数据 提取 
天 体 的 数目 。 从 表 1 中 匹配 的 数目 可 以 看 出 ， 这 种 针对 位 置 的 匹配 度 均 在 85% 以 上 ， 
考虑 到 噪声 影响 和 提取 方法 的 差异 ， 这 种 匹配 程度 是 可 以 接受 的 。 

5 是 对 小 于 29 等 的 天 体 数目 进行 统计 ， 从 该 图 中 可 以 看 到 ， 对 于 模拟 图 像 
所 提取 的 天 体 数 目 和 人 参考 数据 天 体 的 数目 随 着 亮度 变化 的 趋势 是 比较 一 致 的 。 


a. 模 拟 图 像 b. 原始 图 像 
图 4 模拟 图 像 结果 


Fig. 4 Simulation Image vs. Original Image(a: Simulation Image, b: Original Image ) 
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图 5 提取 天 体 数目 统计 对 比 


Fig. 5 Extracted Galaxies from Simulation Image vs. HUDF catalog 


表 1 利用 模拟 图 像 提 取 的 天 体 与 HUDF catalog 天 体 匹 配 数目 统计 


Tab. 1 Extracted galaxies analysisi 


参考 数据 中 天 体 数目 根据 模拟 提取 的 天 体 ”匹配 的 数目 匹配 率 
数据 数目 

<26 mag 匹配 天 体 数 452 406 391 86.5% 
=| 
<27 mag 匹配 天 体 数 871 772 753 86.45% 
<28 mag 匹配 天 体 数 ”1588 1400 1363 85.83% 
目 
<29 mag 匹配 天 体 数 ”2655 2490 2391 90.1% 
目 


图 6 是 将 模拟 数据 提取 出 的 天 体 的 有 效 半径 与 参考 数据 天 体 的 有 效 半径 进行 对 
比 ， 采 用 了 计算 匹配 天 体 有 效 半 径 之 差 的 方法 进行 对 比 ， 匹 配 天 体 的 方法 见 表 1 
说 明 ， 在 这 里 只 是 确定 了 位 置 上 的 匹配 ， 这 两 个 统计 的 样本 都 是 亮度 小 样 29 等 天 
体 。 图 6 (a) 统计 该 差 值 随 着 亮度 分 布 的 变化 ， 横 坐标 是 亮度 ， 纵 坐标 是 参考 数 
据 中 的 天 体 的 有 效 半径 与 相应 的 模拟 数据 提取 的 天 体 的 有 效 半径 之 差 ， 从 该 分 布 
图 上 可 以 看 出 ， 亮 度 越 小 ， 其 有 效 半径 的 差 值 也 越 小 ， 这 一 点 也 符合 我 们 的 认识 。 
图 6 (b) 是 对 差 值 的 分 布 进行 统计 ， 从 这 个 统计 可 以 看 出 ， 这 一 差 值 绝 大 多 数 都 
分 布 在 0 值 左 右 ， 有 部 分 存在 这 较 大 的 偏差 。 导 致 存在 这 种 差异 有 很 多 原因 ， 其 中 
包括 提取 样本 的 精度 (2.1 节 中 对 两 种 参考 样本 的 分 析 ) ， 图 像 模拟 中 加 入 的 信 
号 噪声 、 背 景 噪 声 和 设备 的 噪声 及 利用 SExtractor 对 天 体 进 行 提取 的 参数 也 存在 
一 定 差异 。 
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(a) 有 效 半径 之 差 随 着 天 体 亮 度 的 分 布 (b) 有 效 半径 之 差 的 统计 分 


6 模拟 图 像 所 提取 天 体 的 有 效 半径 分 析 图 


Fig. 6 Anayasis of Effective Radius for Galaxies from Simulation Image 


图 7 是 对 模拟 图 像 所 提取 的 天 体 的 亮度 与 参考 数据 的 天 体 亮度 进行 的 对 比 ， 
对 比 的 方法 同 图 6(a) 一 样 ， 通 过 亮度 之 差 随 着 天 体 亮度 的 分 布 进行 分 析 。 该 图 中 ， 
横 坐 标 是 完 度 ， 纵 坐标 是 参考 数据 中 天 体 的 亮度 与 相应 的 模拟 数据 中 提取 的 天 体 
的 亮度 之 差 。 统 计 的 样本 依然 是 亮度 小 于 29 等 天 体 。 从 该 图 中 可 以 看 出 ， 这 一 亮 
度 差 值 大 部 分 分 布 在 0 轴 的 上 下 ， 偏 差 绝 大 多 数 都 集中 在 0. 5 以 下 。 从 这 张 图 中 ， 
还 可 以 看 出 一 个 规律 ， 对 于 亮度 越 大 的 天 体 ， 模 拟 后 再 提取 ， 其 亮度 值 越 准 确 。 
2.3 性 能 分 析 


本 文中 图 像 模拟 软件 开发 以 及 统计 的 数据 生成 都 是 用 JAVA 程序 开发 。 在 天 文 
中 用 JAVA 开发 的 程序 较 少 ， 本 文 之 所 以 选择 JAVA 程序 开发 主要 只 因为 JAVA 程序 
存在 着 以 下 一 些 优势 : 1 JAVA 去 掉 了 指针 ， 有 自动 的 垃圾 回收 机 制 ， 这 样 大 大 减 
小 编写 程序 时 由 于 内 存 泄露 出 现 的 问题 ， 因 此 ， 能 够 将 更 多 的 时 间 用 在 算法 的 编 
BE. 

2. JAVA 程序 是 面向 对 象 的 编程 ， 在 程序 的 可 读 性 上 要 好 于 C 程序 
3. JAVA 支持 并 行 计 算 ， 且 编写 简单 ， 对 于 本 程序 的 并 行 计算 上 有 很 大 的 帮助 。 

另外 ， 本 文 的 程序 运行 在 Linux 系统 下 ， 硬 件 配置 CPU 为 Intel core i5 4 
核 3. 16Hz， 内 存 为 166B， 在 程序 运行 时 ，4 核 全 部 使 用 ， 计 算 天 体 的 数目 为 
4107， 计 算 所 需要 的 时 间 为 1267s， 平 均 每 个 天 体 所 需要 时 间 为 308 ms。 程 序 运 行 
时 峰值 所 需要 的 内 存 为 4. 16B。 

程序 运行 的 大 部 分 时 间 和 内 存 都 是 消耗 在 利用 Sersic 模型 计算 天 体 的 形态 ， 
在 计算 天 体 的 形态 时 ， 由 于 1. 4 节 中 所 述 子 像素 的 需求 以 及 计算 的 精度 ， 需 要 将 一 
个 像素 又 分 解 成 若干 个 小 像素 ， 计 算 的 时 间 和 内 存 的 消耗 和 像素 分 解 的 个 数 成 正 
比 。 

对 于 星系 形态 的 计算 采用 了 二 维 的 数值 积分 ， 为 了 计算 更 加 准确 ， 本 文 将 一 
ee i me e a 个 大 的 像 元 的 准确 数 
据 ， 文 中 兼顾 性 能 需求 没有 对 像 元 进行 无 限 划分 ， 只 是 分 解 为 32X32 个 子 像素 ， 
4500 X 4500 个 像素 的 图 像 就 变 成 了 144, 000X 144, 000 个 像素 的 图 像 ， 这 样 数据 量 


已 经 很 大 ， 采 用 并 行 及 算法 上 的 优化 使 得 每 个 天 体 的 计算 时 间 达 到 308ms 。 
3 总 结 及 展望 


本 文中 对 于 图 像 模拟 所 用 到 的 方法 都 是 天 文中 比较 基础 的 模型 ， 例 如 Sersic 
模型 、Gaussian 模型 。 将 这 些 基础 模型 应 用 到 模拟 中 ， 可 以 达到 简单 高 效 的 目的 。 
但 是 简单 高 效 只 是 在 模拟 结果 真实 正确 的 基础 上 才能 成 立 。 通 过 本 文 第 3 部 分 的 分 
析 可 以 发 现 ， 模 拟 的 结果 与 真实 数据 存在 一 定 差异 ， 但 是 基本 特征 都 能 够 正确 的 
表现 出 来 〈 例 如 ， 大 小 ， 方 向 、 亮 度 等 ) ; 通过 统计 分 析 可 以 看 出 ， 对 于 模拟 图 
像 提取 出 的 天 体 的 数目 、 亮 度 、 大 小 与 参考 数据 比较 ， 错 误 的 数量 都 在 可 以 控制 
的 范围 内 。 另 外 ， 本 文中 对 模拟 程序 实现 了 并 行 计 算 ， 而 且 在 程序 中 也 做 了 已 做 
了 一 些 优化 《〈 例 如， 三 角 函 数 用 插值 方法 计算 ) ， 通 过 2 部 分 中 对 于 性 能 分 析 可 以 
看 到 ， 计 算 天 体 平 均 时 间 达 到 了 毫秒 级 ， 再 次 说 明了 本 文 的 模拟 成 像 快速 。 

模拟 程序 中 ， 点 扩散 函数 只 用 到 Gaussian 函数 ， 虽 然 简单 ， 但 是 对 于 反映 图 
像 的 真实 性 上 还 存在 着 一 定 差别 。 另 外 ， 模 拟 的 目的 主要 是 为 一 些 新 的 巡天 项 目 
分 析 服 务 。 所 以 本 文 的 下 阶段 目标 就 是 在 模拟 程序 上 对 点 扩散 函数 进行 改进 ， 加 
入 真实 光学 系统 特性 ; 另外 需要 引入 光谱 的 模拟 以 适应 多 波段 成 像 模拟 。 同 时 ， 
在 分 析 上 ， 将 模拟 与 新 的 巡天 项 目 相 结合 ， 得 到 与 不 同 硬件 相关 的 分 析 。 
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图 7 模拟 图 像 中 提取 天 体 亮度 分 析 图 


Fig. 7 Magnitude analysis for Galaxies from Simulation Image 


附录 1. SExtractor 提取 主要 参数 


DETECT_MINAREA 3 # minimum number of pixels above threshold 
DETECT_THRESH 1.2 # <sigmas> or <threshold>,<ZP> in 
mag. arcsec—2 

ANALYSIS THRESH 2 # <sigmas> or <threshold>,<ZP> in 


mag. arcsec—2 


FILTER Y # apply filter for detection (Y or N)? 


FILTER NAME default. conv # name of the file containing the filter 
DEBLEND_NTHRESH 32 # Number of deblending sub-thresholds 
DEBLEND MINCONT 0.005 # Minimum contrast parameter for 
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Abstract: Image simulation has played an increasingly important role in 
astronomical research. There are several important effects as fol low: 
Through simulation, it can provide an important basis for the 
evaluation of astronomical observation equipment; Data processing 
algorithm can be verified through the processing of simulation data. In 
this paper, we describe the whole process of image simulation and 
optimize the local methods of image simulation, so that the simulation 
results can be generated efficiently and accurately. The main work of 
this paper includes the following aspects: Simulation for astronomical 
objects, including point sources (such as stars) and extended source 
(galaxies), and mainly introduce the galaxy simulation in this paper; 
Simulation of observation conditions, mainly include instrument Noise, 
stray light etc. ; Simulation of optical system character which shows in 
image as point spread function (PSF); Optimization of calculation flow 
from calculations and procedural methods; Analysis for extracted 
galaxies sample. The data used in this paper is Hubble Ultra Deep Space 
(HUDF) ACS WFC i-band (F775) data with a limit of magnitude equal to 29 
AB Mag. 
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